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Using an isospin-dependent quantum molecular dynamics, nuclear stopping in 
intermediate heavy ion collisions has been studied. The calculation has been done for 
colliding systems with different neutron-proton ratios in beam energy ranging from 
15MeV/u to 150MeV/u. It is found that, in the energy region from above Fermi 
energy to 150MeV/u, nuclear stopping is very sensitive to the isospin dependence 
of in-medium nucleon-nucleon cross section, but insensitive to symmetry potential. 
From this investigation, we propose that nuclear stopping can be used as a new probe 
to extract the information on the isospin dependence of in-medium nucleon-nucleon 
cross section in intermediate energy heavy ion collisions. 
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Nuclear stopping in heavy ion collisions (HIC) has been studied by means of rapidity distribution!^! 
and asymmetry of nuclcon momentum distribution^^'^'^l. It is an important quantity in determining the 
outcome of a reaction I^'^l. W. Bauer et al. pointed out that in intermediate energy HIC, nuclear stopping 
power is determined by both the mean field and the in- medium nucleon-nucleon (N-N) cross section I^'^l , 
but the mean field he used did not include symmetry potential. Recently S. A. Bass, S. J. Yennello, H. 
Johnston, and Bao-An Li et al suggested that the degree of approaching to isospin equilibration provides a 
means to probe the mechanism and the power of nuclear stopping in HIC But it is not clear how the 

stopping power depends on the symmetry potential in the same collisions. In this Letter, we report a new 
possibility to extract information on the in-medium N-N cross section in intermediate energy HIC by using 
nuclear stopping as a probe. The effects of both in-medium N-N cross section and symmetry potential on 
nuclear stopping shall be studied comparatively for colliding systems with different neutron-proton ratios in 
the beam energy ranging from 15McV/u to 150McV/u by using an isospin dependent quantum molecular 
dynamics (IQMD) model. The calculated results show that stopping power depends strongly on the isospin 
dependence of the in-medium N-N cross section, but weakly on the symmetry potential for all the colliding 
systems studied in the beam energy region from 45MeV/u to 150MeV/u. However, in the energy region 
below Fermi energy, nuclear stopping is sensitive to both the in-medium N-N cross section and the symmetry 
potential. 

The following two quantities can be used to describe nuclear stopping in HIC. The momentum quadrupole 
Qzz defined as Qzz = X^f (2-Pz(i)^ ~ Px{iY ~ -^^(O^)' the transverse-parallel ratio of momentum R 
given by i? = (2/7r)(^^ | P±{i) \)/(^f \ I)- Here the total mass A is the sum of the projectile mass 

Ap and the target mass At. The values of the transverse and parallel components of the momentum of the 
i-th nuclcon are P±{i) = \/ Px{iy + Py{iY and P//{i) = Pzii), respectively. 

In order to describe the isospin effects on the dynamical process of HIC, quantum molecular dynamics 
(QMDI^^I) should be modified properly: (1) the density dependent mean field should contain the correct 
isospin-dependent terms including symmetry potential and Coulomb potential, (2) the in-medium N-N cross 
section should be different for neutron-neutron ( proton-proton ) and neutron-proton collisions, and finally, 
(3) Pauli blocking should be counted by distinguishing neutron and proton. The interaction potential is 
given by 

JJ jjSky _|_ jjYuk _|_ jjCoul _|_ jjMDI _|_ jjPauli _|_ jjSym /-|^\ 
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where U^'^^ is the density-dependent Skyrme potential, 

C/S'=f = a(-^)+/3(-^)T (2) 
Po Po 

a = -390.1Mey, /3 = 2,2Q.^MeV and 7 = 1.1667. C/'^'"'' the Coulomb potential and f/^"*^ the Yukawa 
potential (for detail, see Refs.[18]). JJ^^^ is the momentum-dependent interaction!^ ^1, 

[/^^^ = Uln^MPt - nf + 1]-^ (3) 

Po 

where the parameters t4 = 1.57MeV, = 5 x lQ~^MeV~'^ , p and po are nuclear density and its normal 
value, respectively. J/^"^'* is the Pauli potential, 

UPauH ^ V,(-^)^exp(-M--^ - ^^^^)<5,,p, (4) 
PoQo ^Po 

where 5p^p■ = 1 for neutron-neutron or proton-proton, and 6p^p. = for neutron-proton. The parameters 

Vp = 30MeV, Po = 4:0QMeV/c and qq = 5.64 fm. According to our experience, the interaction potentials 
containing Pauli potential can describe the structure effect of fragmentation in the dynamical process of HIC 

[19] _ 

jjsym jg ^jjg symmetry potential. In the present calculation, three different density dependences of 
the symmetry potential I^^] are used, i.e., U^^"" = cFi{u)St^, C/|^™ = cF2iu)[ST, + ^6^], and C/l"" = 
cFs{u)[5tz — jS"^], where = 1 for neutron and — —1 for proton, c is the strength of symmetry potential, 
taking the value of or 32MeV. -Fi(w) = u, F2{u) = and Fz{u) = v}/"^, u = p/ po- ^ is the relative neutron 
excess 5 = = ^""^^ , where p„ and pp are neutron and proton densities, respectively. First of all, the 

density distribution of nucleus was calculated by using Skyrme-Hatree-Fock with parameter set SKM*P^1. 

An empirical expression of the in- medium N-N cross section is used, (j^'j^ = (1 + ct-^)(j{j^, with 
the parameter a w —0.2, and cr^^ is the experimental free N-N cross scctionl^-'^l. It is known that the free 
neutron-proton cross section is about 3 times larger than the free proton-proton or neutron-neutron cross 
section below 300 MeV. 

The nuclear stopping are studied for the colUding systems ^°iVe -|-^° Ne, *°^r -|-^° Ar, ^^Zn Zn, 
^^^Sn -f"^ Sn and ^^^5n -|-^^^ Sn at different beam energies ranging from 15MeV/u to 150 MeV/u. The 
neutron-proton ratios of the above colliding systems are 1.0, 1.22, 1.67, 1.24, and 1.48, respectively. To make 
a comparative study of the isospin effects of the in-medium N-N cross section and the symmetry potentials, 
we have investigated the following three different cases, (i) The symmetry potential U^y™ = {7*^"* is used 
with the strength c = 32 MeV and the isospin dependent in-medium N-N cross section is employed. This 



case is denoted by C=?7*^™ + cr'*" and solid lines in the figures, (ii) J/*""* = C/^""* and cr]^^ is isospin- 
independent, denoted by C—Ui^"^ + cr""*"" and dash lines, (iii) There is no symmetry potential (c = 0) and 
crj^l^ is isospin-dependent, denoted by C=0 + a^^° and dotted lines. 

Fig.l depicts the time evolution of R for central collisions of the systems ^^'^Sn +^^^ Sn (top row), and 
i^^Sn+i^^Sn (bottom row), with the beam energies of 15MeV/u, 25MeV/u, 35MeV/u, 45MeV/u, 72MeV/u, 
lOOMeV/u, and 150MeV/u for the above three cases. It is very clear that for both colliding systems the 
nuclear stopping R depends strongly on the isospin dependence of the in-medium N-N cross section and 
weakly on the symmetry potential as the beam energy is above 45MeV/u. As the beam energy decreases to 
below the Fermi energy, the nuclear stopping depends on both isospin dependence of in-medium N-N cross 
section and the symmetry potential. 

In Fig. 2, the time evolution of Qzz is shown for central collisions of the system ^^^Sn + ^^^Sn at different 
beam energies in the case of C = f/*^"' + o-*"". It is clearly shown that as decreasing the beam energy, 
the asymptotic value of Qzz decreases towards zero, indicating an isotropic nucleon momentum distribution 
of the whole composite system and consequently a full stopping at the beam energy below Fermi energy. 
As beam energy increases to above Fermi energy, because of the pre-equilibrium particle emission and 
the longitudinal motion of the project-like and target-like nuclei, Qzz will got certain non-vanishing value, 
indicating partial transparency. It is also seen that the relaxation time decreases as increasing beam energy, 
which shows the fact that high beam energy leads to more violent N-N collisions and faster dissipation. This 
is consistent with the isospin equilibrium process as shown by Bao-An Li et al. I^^l . 

In order to trace neutron and/or proton observables, in Fig. 3 is given the time evolution of Qzz for 
neutron and proton respectively. The conclusion drawn from Fig. 1 is applicable to both neutron and proton 
stoppings. 

In Fig. 4 is plotted the time evolution of R for the reaction ^^^Sn -|- ^^^Sn at E = lOOMeV/u for different 
impact parameters. It is noticed that the nuclear stopping for small impact parameters (b=O.Ofm, l.Ofm, 
2.0fm) shows the same behavior as the central collisions as in Fig. 1, while as increasing the impact parameter, 
the behavior becomes different and the dominant role played by the isospin dependent in-medium N-N cross 
section gradually disappears. The impact parameter dependence of R{b) can be seen more clearly from Table 
1 where the asymptotic values of R are given for different b. 

Table 1 The asymptotic values of the stopping R for the reaction ^^"^ Sn+^^'^ Sn 
at E=100MeV/u as function of impact parameter 



4 



b(fm) 


0.0 


1.0 


2.0 


3.0 


4.0 


5.0 




0.70 


0.72 


0.60 


0.55 


0.51 


0.50 




0.82 


0.82 


0.71 


0.63 


0.56 


0.50 


i?(0 + cr'"°) 


0.84 


0.81 


0.74 


0.64 


0.58 


0.51 



In Fig. 5 is shown the impact parameter-averaged asymptotic vahics ( the values at t > 200fm/c when 
the nuclear stopping becomes nearly a constant as shown in Figs. 1 and 2) of Qzz per nucleon (left window) 
and R (right window) as a function of At + Ap in the following seven cases: C=0 + a"" by the solid square; 
Q^jjsym _^ ^iso denoted by the solid circles; G=U^^"' + a''" by solid triangle; G=U^^"' + a''" by solid 
diamond; 0=01^"^ + cr"""" by open circles; C=f/2^'" + <j"°**'° by open triangle; C=[/3^'" + 0-"°*'*° by open 
diamond. In the calculation the beam energy is 100 MeV/u and five colliding systems have been considered, 
i.e., ^OATe +20 Ne, ^°Ar +4° Ar, ^"^Zn Zn, ^^"^Sn +"2 Sn and ^^^Sn +^^^ Sn. From Fig. 5, it can be seen 
that the nuclear stopping depends strongly on the isospin dependence of in-medium N-N cross section and 
weakly on the symmetry potential, though it is slightly different for the three different forms of symmetry 
potential, especially for larger-mass systems. It is also shown that the nuclear stopping power increases as 
increasing the total mass of the colliding system, the decrease of Qzz or the increase of R implies the increase 
of nuclear stopping. 

Fig. 6 shows the asymptotic value at t = 200fm/c of i? as a function of the beam energy for the same 
colliding systems and the same collision conditions as in Fig. 1. From the figure, one can draw the same 
conclusion as from Fig. 1. 

From the above comparative study of nuclear stopping for different colliding systems, beam energies, and 
impact parameters by using IQMD, we have obtained the following physical picture of the collision dynamics: 
from above the Fermi energy to about 150MeV/u, the dynamics is dominated by N-N collisions and the role 
played by the mean field is less important. The main consequence of the N-N collisions is the transformation 
of the initial longitudinal motion to the motion in other directions and the subsequent thermalization of 
the system. In this case, nuclear stopping R and Qzz are sensitive to nucleon- nucleon collisions, and can 
thus be used as a measure of the dissipative process. However, as the beam energy decreases to below the 
Fermi energy, the collision dynamics is governed by both the mean field and the N-N collision which result in 
one-body and two-body dissipation or thermalization, respectively. Thus as indicators of the thermalization, 
R and Qzz depend on both the mean field and the in-medium N-N cross-section. 

In summary, from the above results and discussions we conclude that the nuclear stopping R and Qzz 
can be used as a probe to extract information on the isospin dependence of the in-medium N-N cross section 
in HIC in the beam energy region from above the Fermi energy to about 150 MeV/u. 
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Figure 3: The time evolution of momentum quadrupoles for neutron Qzz{neutron) and proton Qzz{proton) 
for the central coUisions of ^^'^Sn +^^^ Sn at the beam energies of 35MeV/u, 72MeV/u, and lOOMeV/u in 
the three cases as shown in Fig.l. 
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